Particle-Hole Optical Model: Fantasy or Reality? 
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An attempt to formulate the optical model of particle-hole-type excitations (including giant res- 
onances) is undertaken. The model is based on the Bethe-Goldstone equation for the particle-hole 
Green function. This equation involves a specific energy-dependent particle-hole interaction that 
is due to virtual excitation of many-quasiparticle configurations and responsible for the spreading 
effect. After energy averaging, this interaction involves an imaginary part. The analogy between 
the single-quasiparticle and particle-hole optical models is outlined. 
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Damping of giant resonances (GRs) is a long-standing 
problem for theoretical studies. There are three main 
modes of GR relaxation: (i) particle-hole (p-h) strength 
distribution (Landau damping) , which is a result of the 
shell structure of nuclei; (ii) coupling of (p-h)-type states 
with the single-particle (s.p.) continuum, which leads to 
direct nucleon decay of GRs; and, (iii) coupling of (p- 
h)-type states with many-quasiparticle configurations, 
which leads to the spreading effect. An interplay of these 
relaxation modes takes place in the GR phenomenon. 

A description of the giant-resonance strength function 
with exact allowance for the Landau damping and s.p. 
continuum can be obtained within the continuum-RPA 
(cRPA), provided the nuclear mean field and p-h inter- 
action are fixed [l|. As for the spreading effect, it is 
attempted to be described together with other GR re- 
laxation modes within microscopic and semimicroscopic 
approaches. The coupling of the (p-h)-type states, which 
are the doorway states (DWS) for the spreading effect, 
with a limited number of 2p-2h configurations is explic- 
itly taken into account within the microscopic approaches 
(see, e.g., Refs. [2|,[3j]). Some questions to the basic points 
of these approaches could be brought up: (i) "Thermal- 
ization" of the DWS, which form a given GR, i.e. the 
DWS coupling with many-quasiparticle states (MQPS) 
(the latter are complicated superpositions of 2p-2h, 3p- 
3h, ...configurations), is not taken into account. As a 
result, each DWS may interact with others via 2p-2h 
configurations. Due to complexity of MQPS one can rea- 
sonably expect that after energy averaging the interac- 
tion of different DWS via MQPS would be close to zero 
(the statistical assumption), (ii) The use of a limited 
basis of 2p-2h configurations does not allow to describe 
correctly the GR energy shift due to the spreading effect. 
The full basis of these configurations should be formally 
used for this purpose, (iii) With the single exception of 
Ref. [4j , there are no studies of GR direct-decay proper- 
ties within the microscopic approaches. 

Within the so-called semimicroscopic approach, the 
spreading effect is phenomenologically taken into account 
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directly in the cRPA equations in terms of the imaginary 
part of an effective s.p. optical-model potential [9, Q- 
Within this approach, the afore-mentioned statistical as- 
sumption is supposed to be valid and used in formulation 
of the approach. The GR energy shift due to the spread- 
ing effect is evaluated by means of the proper dispersive 
relationship, and therefore, the full basis of MQPS is for- 
mally taken into account [7j . The approach is applied to 
description of direct-decay properties of various GRs (the 
references are given in Ref. [6|). In accordance with the 
"pole" approximation used for description of the spread- 
ing effect within the semimicroscopic approach, the latter 
is valid only in the vicinity of the GR energy. However, 
for an analysis of some phenomena it is necessary to de- 
scribe the low- and/or high-energy tails of various GRs. 
For instance, the asymmetry (relative to 90°) of the (7n)- 
reaction differential cross section at the energy of the 
isovector giant quadrupole resonance is determined, in 
particular, by the high-energy tail of the isovector giant 
dipole resonance |8|. Another example is the isospin- 
selfconsistent description of the IAR damping. In partic- 
ular, the IAR total width is determined by the low-energy 
tail of the charge-exchange giant monopole resonance [9j- 

El. 

In the present work, we attempt to formulate a model 
for phenomenological description of the spreading effect 
on p-h strength functions at arbitrary (but high enough) 
excitation energies. The formulation of this semimicro- 
scopic model (simply called as the p-h optical model) 
is analogous to that of the single-quasiparticle optical 
model [la . ll3j . The dispersive version of this model 12] is 
widely used for description of various properties of single- 
quasiparticle excitations at relatively high energies (see, 
e.g., Ref. [ll|). This model can be also used for descrip- 
tion of direct particle decay of subbarrier s.p. states [15j . 

The starting point in formulation of the single- 
quasiparticle optical model is the Fourier-component 
of the Fermi-system single-quasiparticle Green function, 
G(x,x';e), taken in the coordinate representation (see, 
e.g., Refs. [12, LL3|])- By analogy with that, we start 
formulation of the p-h optical model from the Fourier- 
component of the Fermi-system (generally, nonlocal) p- 
h Green function, A(x,x';xi,x' 1 ;lo), also taken in the 



coordinate representation. Being a kind of the Fermi- 



system two-particle Green function (definitions see, e.g., 
in Ref. 116]), A satisfies the following spectral expansion: 



A(x,x'; xi,x'i',u) = 
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Here, uj s = E s — Eq is the excitation energy of an exact 
state \s) of the system and p s (x,x') = (s\^ + (x)^>(x')\0) 
is the transition matrix density (^ + (x) is the operator of 
particle creation at the point x). In accordance with 
the expansion of Eq. ([1]) the p-h Green function de- 
termines the strength function Sy(ui) corresponding to 
an external (generally, nonlocal) single-quasiparticle field 

V = J $+(x)V(x, x')${x')dxdx' = [$+^§ 



and 



S v (u) = lm\V+A(w)V] , 

n L J 



(2) 



where the brackets [. . . ] mean the proper integrations. 

The free s.p. and p-h Green functions, Gq(x, x'\e) 
and Ao(x, x';Xi,x'i;oj), respectively, are determined by 
the mean field (via the single-quasiparticle wave func- 
tions) and the occupation numbers (only nuclei without 
nucleon pairing are considered). Being determined by 
Eq. (UJ, the free transition matrix densities p s (x, x') are 
orthogonal: p s p , — 5 S8 '. As for the transition den- 
sities p s (x = x'), which appear in the spectral expan- 
sion for the free local p-h Green function Aq(x, Xi] of) — 
Ao(x = x',Xi = x'i,ijS), this statement is wrong. The 
RPA p-h Green function, Arpa(x, x';x±, x'^oj), is deter- 
mined also by a p-h (local) interaction F(x, x'; xi,x[) — 
F(x, x\)5(x — x')5{x\ — x'i), which is responsible for long- 
range correlations leading to formation of GRs. In partic- 
ular, the Landau-Migdal forces F(x, Xi) — > F(x)5(x—x 1 ) 
are used in realizations of the semimicroscopic approach 
of Refs. 0, Q. The RPA p-h Green function satis- 
fies the expansion, which is similar to that of Eq. ([TJ. 
In such a case, the RPA states \d) are the DWS for 
the spreading effect. The local RPA p-h Green func- 
tion Arpa (x, Xi ; co) determined by the p-h interaction 
F(x, X\) is used for cRPA-based description of the GR 
strength function corresponding to a local external field 
V = J y + (x)V(x)$(x)dx Q. 

The s.p. and p-h Green functions satisfy, respectively, 
the Dyson and Bethe-Goldstone integral equations: 



A(u) = A rpa {uj) + [Arpa(w)V(u)A(w)] , (4) 



where 

A RPA (u) = A q (lu) + [Ao(w)TArpa(w)] ■ (5) 
The self-energy operator T,(x,x';e) and the specific p- 
h interaction (polarization operator) V(x, x';xi,x'i, oj) 
describe the coupling, correspondingly, of single- 
quasiparticle and (p-h)-type states with proper MQPS. 
Analytical properties of G and S are nearly the same. A 
similar statement can be made for A and V '. The quan- 
tities S(e) and V(uj) both exhibit a sharp energy depen- 
dence due to a high density of poles corresponding to 
virtual excitation of MQPS. Concluding consideration of 
the basic relationships given above in a rather schematic 
form, we present the alternative equation for the p-h 
Green function: 



A(u) = Ao(u) + [A)M (J 7 + V{u)) A(u)} , 



(6) 



G(e) = G (e) + [G (e)E(e)G(e)] 



(3) 



which follows from Eqs. (@}, ([5]). 

Since the density of MQPS, p m , is large and described 
by statistical formulae, only the quantities S(x, x'\ e) and 
P(x,x';Xi,x'i;oj) averaged over an interval J 3> p^ 1 
can be reasonably parameterized. As applied to E(e) = 
S (e + iJSgn(e — p,)), it is done, e.g., in Refs. (l2i.ll3j: 

S(x, x';e) — Sgn(e — p) {—iw(x; e) + p(x; e)} 6(x ~ x'). 

(7) 
Here, p is the chemical potential and w(x; e) is the imagi- 
nary part of a (local) optical-model potential. Assuming 
that the radial dependencies of p and w are the same, 
i.e. w(x;e) — > w(r)w(e) and p(x;e) —> w(r)p(e), the in- 
tensity of the real addition to the mean field, p{e), has 
been expressed in terms of w(e) via the corresponding 
dispersive relationship [lj]. It is noteworthy, that the 
optical-model addition to the mean field can be taken as 
the local one, i.e. S(x,x';e) ~ 5(x—x'), in view of a large 
momentum transfer (of order of the Fermi momentum) 
at the "decay" of single-quasiparticle states into MQPS. 
The energy averaged single-quasiparticle Green function 
G(x, x';s) satisfies the Eq. (O, which involves in such a 
case the quantity E of Eq. (J7]) . Actually, G is the Green 
function of the Schrodinger equation, which involves the 
addition to the mean field considcratcd above. 

The energy-averaged polarization operator can be pa- 
rameterized similarly to Eq. ([7]): 



P(x,x';xi,x' 1 ;uj) = {—i\V(x,x';uj) + P(x,x';w)}S(x — Xi)S(x' — x[). 



(8) 



Assuming that the coordinate dependencies of the quan- 
tities P and W are the same, i.e. W(x, x';oj) — > 
W(x, x')W(uS) and P(x,x';oj) -)■ W(x,x')P(u), we can 
express P(u>) in terms of W(u>) via the corresponding dis- 
persive relationship. The example of such a relationship 
is given in Ref. 0. In accordance with Eqs. flU), ©, 
((HJ the energy-averaged p-h Green function satisfies the 
equivalent equations: 

A(u) = A rpa (oj) + [A R pa{u)V{oj)A{u))] , (9) 



A{w) = A>H + [Ao(u) {F + ?(«)) ■*(«)] ■ ( 10 ) 



Formally, Eqs. (j8|)- (fT0| are the basic equations of the 
p-h optical model. In particular, the energy-averaged 
strength function is determined by Eq. ([2]) with the sub- 
stitution A(uj) — > A(u>). 

To realize the model in practice, a reasonable 
parametrization of Imp should be done with taking the 
statistical assumption into account. For this purpose, 
we consider the quantity Arpa within the discrete-RPA 
(dRPA) in the "pole" approximation. In accordance with 
Eq. (TTJ), we have 



A 



RPA 



(x,x';xi,x[;u>) 



d 



P*d( x '' X )Pd(xi,Xl) 

ui — ujd + iO 



(11) 



The statistical assumption [p d Vpd'] ~ S dd > is fulfilled, 
provided that: (i) the intensity W(x, x';u>) is nearly con- 
stant within the nuclear volume, i.e. W(x, x';u>) — > 
W(w); and, (ii) the dRPA transition matrix densities are 
orthogonal, i.e. [p d Pd'} — $dd'- Under these assumptions, 
the solution of Eq. © can be easily obtained in the pole 
approximation: A(ui) = Arpa{w + iW(u>) — P (w))- As 
a result, the energy-averaged strength function is the su- 
perimposition of the DWS resonances: 



Sv(oj) 



1 



-Im 



[Vpd] 



^ cu - cu d + iW{uj) - P{uj) ' 



(12) 



The quantity 2W can be considered as the mean DWS 
spreading width (r^), which might be larger than the 
mean energy interval between neighboring DWS reso- 
nances. 

A few points are noteworthy in conclusion of the above- 
given description of the p-h optical model. Within the 
model the spreading effect on formation of (p-h)-type 
excitations is described phenomenologically in terms of 
the specific (p-h) interaction V. Because the interfer- 
ence between the spreading of particles and holes is taken 
into account by this interaction, the latter cannot be 
expressed via the single-quasiparticle self-energy opera- 
tor 2. Formally, the p-h optical model is valid at ar- 
bitrary (but high enough) excitation energy. The low 



r 



limit is determined by the possibility of using the sta- 
tistical formulae to describe the MQPS density. Within 
the semimicroscopic approach, the substitution like uj — > 
cj + iW(w) — P(w) is used in the cRPA equations to take 
the spreading effect phenomenologically into account in 
the "pole" approximation together with the statistical 
assumption [6|, 0|- Thus, the parametrization of W(w) 
can be taken in the form widely used for the intensity 
of the imaginary part of the effective s.p. optical-model 
potential in implementations of the semimicroscopic ap- 
proach. Within the s.p. optical model the statistical 
assumption for "decay" of different single-quasiparticle 
states with the same angular momentum and parity into 
MQPS seems to be valid. At high excitation energies 
|e — /i |, when the empirical value of w(s) is compara- 
ble with the energy interval between the afore-mentioned 
single-quasiparticle states, the empirical radial depen- 
dence w(r) becomes nearly constant within the nuclear 
volume (see, e.g., Refs. [14]). 

The p-h optical model can be simply realized in 
terms of the energy-averaged local p-h Green function 
A(x, x\ ; oS) to describe the strength function of a "single- 
level" GR, because in such a case there is no need for the 
statistical assumption. Being more simple, the equations 
like flU-©, ©, (JTOl) are actually the straight-forward 
extension of the corresponding cRPA equations. In prac- 
tice, within the cRPA it is more convenient to use the 
equation for the effective field V{x, uj), which corresponds 
to a local external field V{x) and is determined in accor- 
dance with the relationship: [FA(w)] = [F(w)j4o(w)]- 
The effective field determines the strength function: 



S v {u) = Im [VoM">)V{uj\ , (13) 

and satisfies the equation: 

V(uj) = V+[(F + n(uj))A Q (uj)V(uj)]. (14) 

The energy-averaged local polarization operator can be 
parameterized similarly to Eq. ([%]): 

3(a;,a;i;w) = C{-iW(x;u) + P(x; u)}S(x - xi). (15) 

Here, C — 300MeVfm 3 is the value often used in 
parametrization of the Landau-Migdal forces; W and P 
are the dimensionless quantities, which can be parame- 
terized as follows: W(x;uj) — > W(r)W(uj) and P(x;tu) — > 
W(r)P(oj), where P(ui) is determined by W(u>) via the 
corresponding dispersive relationship [7|. 

Due to strong coupling with s.p. continuum, the high- 
energy GRs (they are mostly the overtones of correspond- 
ing low-energy GRs) can be roughly considered as the 
"one-level" ones. Being the IAR overtone, the charge- 
exchange (in the /^"-channel) giant monopole resonance 
(GMR(') is related to these GRs. Within the isospin- 
selfconsistent description of the IAR damping [10|, lllj . 




Figure 1. The "Coulomb" strength function calculated for 
the 208 Pb parent nucleus within the cRPA (thin line), the 
semimicroscopic approach (dash-dotted line), and p-h the 
optical model (full line). The arrow indicates the IAR energy. 



the low-energy "tail" of the GMR/ ' in the energy de- 
pendence of the "Coulomb" strength function S^ (oj) 
determines the IAR total width Ta via the nonlinear 
equation: 



-io(-), 



T A = 2vS A :'S y c > {Lu = Lo A ) 



(16) 



Here, Sa — {N — Z) is the IAR Fermi strength, uja 
is the IAR energy, and the "Coulomb" strength func- 
tion corresponds to the external field V(x) — > Vq = 
(Uc(r) —ua + %Fa)t(~>, where Uc(r) is the mean 
Coulomb field. Strength function Sq (oj) exhibits a wide 
resonance corresponding to the GMR(~). In Fig. [TJ we 
present the strength function calculated for the 208 Pb 
parent nucleus within: (i) the cRPA (in such a case the 

strength function S@ (u) = uja) determines the IAR 



total escape width found without takin g t he isospin- 
forbidden spreading effect into account J17j|): (ii) the 
semimicroscopic approach [Tl| and, (iii) the p-h opti- 
cal model by Eqs. (jT5|) — (fT5 |) . All the model parameters, 
parameterization of the imaginary part of the effective 
single-quasiparticle optical-model potential I(r; ui) jllT ] 
and parameterization of W(r;uj) in Eq. (TT51) are taken 
the same in both approaches. The intensities of I(r; to) 
and W(r; to) are chosen to reproduce in calculations the 



observable total width of the GMR' 



*Bi (~ 15 



MeV). Both approaches lead to the similar results, which 
are not exactly the same for the low-energy "tail" of the 
GMR( - ) at u> ~ uja- Irregularities in the energy depen- 
dence of Sq'(uj) calculated within the p-h optical model 
are explained by the fact that the GMR(~) can be roughly 
considered as the "one-level" one. 

In the present work, the optical model of particle- 
hole-type excitations has been formulated in terms of 
the energy-averaged nonlocal particle-hole Green func- 
tion. The equation for this Green function involves a spe- 
cific energy-dependent particle-hole interaction, which is 
due to virtual excitation of many-quasiparticle configu- 
rations. The intensity of the imaginary part of this in- 
teraction should be taken nearly constant within nuclear 
volume to satisfy the statistical assumption on the inde- 
pendent spreading of different particle-hole-type states 
which form a given giant resonance. The strength func- 
tion of the "single- level" giant resonance can be described 
in terms of the energy-averaged local particle-hole Green 
function. 

Along with numerical realizations, the particle-hole op- 
tical model can be extended to describe direct particle 
decays of giant resonances. These points are under con- 
sideration. 

In conclusion, one can say that in formulation of the 
particle-hole optical model we are on a way from fantasy 
to reality. 
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